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We have studied a model of a complex fluid consisting of particles interacting through a hard 
core and a short range attractive potential of both Yukawa and square-well form. Using a hybrid 
method, including a self-consistent and quite accurate approximation for the liquid integral equation 
in the case of the Yukawa fluid, perturbation theory to evaluate the crystal free energies, and mode- 
coupling theory of the glass transition, we determine both the equilibrium phase diagram of the 
system and the lines of equilibrium between the supercooled fluid and the glass phases. For these 
potentials, we study the phase diagrams for different values of the potential range, the ratio of the 
range of the interaction to the diameter of the repulsive core being the main control parameter. Our 
arguments are relevant to a variety of systems, from dense colloidal systems with depletion forces, 
through particle gels, nano-particle aggregation, and globular protein crystallization. 
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I. INTRODUCTION 

Recently novel results in the statistical mechanics of 
the fluid state have emerged from a series of studies of 
the phase diagram of particles with a hard core and short- 
ranged attractive potential. At first sight this type of in- 
teraction potential would have been expected to yield the 
well known phase diagrams of simple solids, liquids and 
gases. Here we shall see that a new scenario emerges. 
In many practical situations attractive forces that are 
short-ranged compared to the size of the particles arise 
because the size of the particles itself is large, whilst 
the physical forces retain their typical microscopic range. 
Thus, in order to model large molecules, such as proteins, 
colloids and nano-particles, one often works in a regime 
where the ratio of the range of attraction to the size of 
the repulsive core is small. This crucial issue is emerging 
in the literature in such areas as protein crystallization 
dense colloids |^-|^, nano-particle assemblies, pre- 
ceramic particle gelation, latex formation, "buckyballs" 
and many others. In some arenas its relevance begins 
to be recognized, while in others it still remains to be 
understood. 

Atoms or molecules interacting via a potential with a re- 



pulsive core and an attractive tail may form gases, liquids 
and solids as a function of the temperature and density. 
The attraction is responsible for the liquid-gas transi- 
tion, while the solid is dominated by the repulsion. In 
the low-to- medium range of densities, particle entropy 
and energy enter the free energy and compete to pro- 
duce liquid (energy-favoured) and gas (entropy-favoured) 
phases. Thus, beneath the critical point, the total sys- 
tem free energy is optimized if the system splits into two 
sub-systems, and thereby a sum of two independent free- 
energies. One of these free energy terms is dominated 
by the energy (liquid) and the other (gas) by the en- 
tropy. A homogeneous, unseparated system, maintained 
at the average density of the two sub-systems, would not 
be able to fully optimize either energy or entropy. This, 
in essence, is the cause of typical first-order phase transi- 
tions between a liquid and a gas. In the liquid, as the den- 
sity increases the remaining diffusive motions of the par- 
ticles, are reduced and their attendant (configurational) 
entropy diminishes. Now, the random structure of a liq- 
uid is favored provided overall diffusive particle motion is 
possible. When, at higher densities, diffusion is greatly 
limited, the free volume is better utilized by making a 
regular crystalline array, and the entropy of the system is 
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increased by the vibrations of particles within the regular 
array of a crystal. This is the reason why hard-sphere flu- 
ids begin to form a regular crystal at about 49% volume 
fraction, and above the freezing transition at 55% have 
a higher entropy than a liquid-like structure at the same 
volume fraction. Indeed, for hard spheres, at volume 
fraction between 49% and 55% the system again phase 
separates to achieve the optimal value of entropy: a lower 
density fluid of 49% volume fraction and a well-packed 
crystal at 55%. Nevertheless, it is well known that in 
some experimental studies the system may be unable to 
access the crystalline state rapidly enough, and the dis- 
ordered super-cooled liquid structure freezes into a glass. 
For hard sphere fluids, this occurs at a volume fraction 
higher than about 58%. The formation of the glass is a 
signal that equilibrium quantities such as the free energy 
may not be sufficient, to describe the behavior of the sys- 
tem, and that long-lived dynamically arrested states may 
be important under such conditions. 

It is possible to relate these qualitative comments to 
recent understanding acquired from experiments, simu- 
lations and theory. Thus, there is a separation of time 
scales in dense systems, including super-cooled liquids, 
in which particles spend a long time trapped by a sur- 
rounding cage and thereafter they escape from it. In fact, 
dynamically slowed systems begin to exhibit a plateau in 
the self-correlation function, reflecting the time spent in 
cages, and then a decay (a-relaxation) reflecting escape 
from the cage and free configurational motion. This ap- 
proach to the glass transition was formalized ^ with 
the idea that motions in dense fluids can be divided into 
intra- and inter-basin motions, where the basins refer to 
the multidimensional potential energy as a function of the 
particles coordinates. This separation of time scales, and 
thereby type of motion, is reasonable for dense systems, 
and may be used to justify a conceptual partitioning of 
the entropy into two parts, a configurational and a local 
contribution |^ . These ideas can be re-expressed by con- 
sidering the system, at a fixed average density, composed 
of central particles trying to escape their cages of neigh- 
bours, which are themselves fluctuating, and exchang- 
ing with their neighbours. Evidently, the a-relaxation 
process corresponds, in systems very close to arrest, to 
the escape of the central particles from their cages af- 
ter some time. In the absence of attractive forces, or 
for relatively long-ranged ones, the motions of the cage 
are restricted by the packing forces. This picture is rel- 
atively clear, at least in a phenomenological manner, for 
hard core spherical particles. For systems that possess a 
strong repulsion, and a long-ranged attraction, there are 
no new special features. Thus both "central" and "cage" 
particles remain in their mutual range of attraction while 
the structural rearrangement takes place. The resulting 
cage breaking can be viewed as being almost the same 
as that for hard spheres, but with changed zero of en- 
ergy. Even if we acknowledge that there are effects due 



to attractions, they can still be considered to be weak 
perturbations of the picture arising from the hard core 
spherical particles. 

Now let us turn from the scenario just described, typ- 
ical of hard-sphere particles or particles where the at- 
tractions are long-ranged compared to the core size, to 
situations where attractions play a principal role. In this 
case, the freedom of the cage particles is considerably 
reduced. Indeed, particles must remain within a cer- 
tain distance from each other where substantial attrac- 
tive energies are still available, or they lose the advan- 
tages of being in the liquid-like structure. This being so, 
the cage around the central particle is much more rigid, 
since only smaller excursions from average positions of 
these cage particles are possible. In fact, for sufficiently 
narrow wells, the time that a particle spends inside the 
cage increases; a plateau regime results, and the system 
will eventually freeze. This "cage rigidity" was also the 
determining factor in hard core particles at very high 
density, leading to the typical colloidal "repulsive glass" . 
Here, though the mechanism and detailed laws will be 
different, we see that short-ranged attractions are able 
to cause cage rigidity, and formation of a solid, either 
a glass, or a crystal. Following these arguments we see 
that for particles with short-ranged attractions, the den- 
sity and temperature window over which the liquid state 
is stable is greatly reduced since the configurational en- 
tropy is reduced. The system separates into a low-density 
fluid state where there is sufficiently free motion, and a 
dense state where there is in any case not sufficient con- 
figurational entropy to sustain a liquid, and the system 
freezes. This state will be a crystal or a dynamically ar- 
rested state. Based on the developing views of others, 
and what we shall present here, all these expectations 
are borne out from precise calculations, and increasingly 
in experiments. 

From the discussion above, we may expect that when 
attractions are short-ranged the arrested glass is of dif- 
ferent nature than the typical repulsive one, since it is 
favoured by both the energy (because the particles are 
not close enough to sample mainly repulsive energy), and 
the local entropy (given that the density is yet relatively 
low) . We have earlier called this arrested state an attrac- 
tive glass and the crystal, previously found by a 
number of others [p^[-p^, we name the attractive crys- 
tal. This distinguishes it from the typical face centred 
cubic (FCC) crystal formed by repulsive forces and the 
analogous repulsive glassy state. 

Given the possibility that for short-ranged attrac- 
tions, solids, both crystals and glasses, can be formed by 
these two distinct mechanisms leading to cage-rigidity, 
we may suppose that it is possible in principle for them 
to co-exist. For fixed (short) ranged potential one way 
of changing the balance of attraction and repulsion is 
to change the density, and we comment that such co- 
existence has been previously shown for crystals [p^|-p^ , 
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and glasses |l^JTl|], and is reproduced here. 

In summation, we expect that, as the well width is 
narrowed, the liquid state becomes progressively less 
favoured, and is replaced by a conventional repulsive 
crystal, or its equivalent glass. An attractive crystal, and 
its equivalent repulsive glass, should also be present at 
higher densities, both types of solid co-existing at some 
typical densities where attractions and repulsions com- 
pete. From the preceding considerations we expect that 
these predictions should be general, irrespective of the 
detailed shape of the potential, and reflecting the typical 
range of the potential. 

We have studied two potential energy models typically 
used to mimic colloidal interactions, the square well (SW) 
and the hard-core Yukawa potential, using a variety of 
techniques of condensed-matter theory. We have deter- 
mined the phase diagrams calculating the liquid free en- 
ergy using perturbation theory for the SW and the self- 
consistent Ornstein-Zernikc approximation (SCOZA) for 
the Yukawa potential. The crystalline free energy has 
been calculated applying second order perturbation the- 
ory for both potentials. This technique has been applied 
to short-ranged potentials before [^-16|, with quite re- 



markable success in reproducing the crystalline free en- 
ergy in comparison to Monte Carlo experiments. How- 
ever, it has been noted that the same method is not 
so satisfactory in calculating the free energy of the liq- 
uid and gas states for very narrow attractive ranges, and 
this affects the accuracy of the phase diagram ]l^ , p7t . 
To avoid this problem in the case of the Yukawa poten- 
tial, we have modified the calculation so that the proven 
good features of the perturbation theory of the crystal 
are combined with the proven quality of the SCOZA for 
the liquid and gas free energies to produce a phase dia- 
gram that is of uniformly good quality. Currently there 
is no working SCOZA method for the SW case, though 
there is no fundamental barrier to developing one. 

In our work the ideal Mode Coupling Theory (MCT) 
Ipsi of super-cooled liquids has been used in order to lo- 
cate the glass transition curves. For colloids this method 
has also been found to describe many elements of the 
transition to the arrested state [|l9[|2C| ], though it does 
suffer from some limitations, a matter to which we return 
later. In the MCT calculations we have used the struc- 
ture factors from Percus-Yevick approximation (PYA) for 
the case of SW and from SCOZA for the case of the 
Yukawa potential. These have proved to be quite accu- 
rate theories of the liquid and fluid states. The SCOZA 
method in particular has been compared to Monte Carlo 
simulations for a range of screening parameters ("well- 
widths" ) of the Yukawa potential, and it has been shown 
pl| , p2t that the agreement for the phase diagrams is 
quantitative, at least for modest values of screening pa- 
rameter. 

By combining the results from the different methods 
described above, we are able to give, for some regimes. 



what we believe to be quite accurate phase diagrams. For 
the general case we believe that the results are at least 
qualitatively correct, and provide us a coherent picture of 
the connection between well- width and the arrangement 
of gas, liquid, crystal and glass phases. It is this overview 
of how the various phenomena fit together that is cur- 
rently missing, and that should prove useful in the vari- 
ous applications alluded to above. We also point out that, 
from experiments, it has become clear that the interac- 
tion between globular proteins, in the range where they 
may crystallize, is characterized by short range attrac- 
tions In the case of proteins, a better comprehension 
of the phenomena would imply the study of anisotropic 
types of potential [^, due to hydrophilic-hydrophobic 
patching of the protein surface. The usual approach is, 
however, to use an effective isotropic interaction obtained 
by averaging over the anisotropy. Thus, as discussed 
later, most of our conclusions will be also important to 
understand the process of protein crystallization. 

The paper is organized as follows. In the next section 
H we describe the approximate closures to the Ornstein- 
Zernike liquid integral equation that we used for the 
Yukawa and SW potentials. The methods employed to 
determine the equilibrium fluid and solid phases are de- 
scribed in sections [II and IV . Section ^ is devoted to 
a brief sketch of MCT applied to attractive potentials, 
while in section VI the relevance of the spinodal curve 
for colloidal systems is discussed. The complete phase 
diagram, including the structural arrest lines is discussed 
in section 



VII, while section VIIl is devoted to our con- 



clusions. 



II. APPROXIMATIONS TO THE OZ EQUATION 
FOR THE YUKAWA AND SQUARE WELL 
MODELS 



In this section we shall discuss the theory used in the 
investigation of the phase diagram. The Ornstein-Zcrnike 
(OZ) equation for the pair correlation function h{r) is 



h{r) = c(r) + p dr c(|r — r \)h{\r |) 



(1) 



where g(r) = h{r) + 1 is the radial distribution function 
and c(r) the direct correlation function. Another impor- 
tant quantity is the static structure factor Sq, which is 
the equal time correlation function of the density vari- 
ables in wave vector space. 



Sq = {p^q{t)pq{t))/N 



(2) 



where the average (• • •) is performed at equilibrium and 
the density variables are Pq{t) = where the 

sum runs over all N particles in the system. The Fourier 
transform of the correlation function hq is related to this 
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quantity by the relation Sq = 1 + phq. The OZ relation 
in the wave vector space reads 



Sa = 



1 



1 - pCq 



(3) 



Cq being the Fourier transform of the direct correlation 
function. 

As it stands Eq. is not closed and some type of 
approximation is needed in order to solve it. In what fol- 
lows we have chosen to calculate structural and thcrmo- 
dynamical properties using the PYA for the SW poten- 
tial and the SCOZA [|l]j2§|2|l for the Yukawa potential. 
Both these model potentials possess some fundamental 
properties that make them good candidates for study- 
ing properties of attractive colloidal particles when the 
range of interaction is short. The use of SCOZA has 
been justified by the success of such approach in predict- 
ing simulation data. 



A. The SCOZA for the Yukawa potential 

The application of SCOZA to a hard-core Yukawa fluid 
provides a semi-analytic calculation of the thermody- 
namic properties of the fluid, liquid and gas states of 
the system pl| , p2t . SCOZA has been applied to Yukawa 
systems with relatively large values of the range of the 
potential, with a satisfactory reproduction of the liquid- 
vapour binodal curves and a good description of the crit- 
ical point region. 

The hard-core Yukawa fluid is described by the follow- 
ing inter-particle potential, 



v{r) 



oo r < a 



ae ■ 



r > a 



(4) 



The parameter e defines the energy scale, while the 
parameter 6, known as screening parameter, determines 
the range of the potential. The larger the b, the shorter 
is the range of the potential. In this paper we set a = 1 
and e = 1, therefore the screening parameter is in units 
of the reciprocal of the hard-core diameter, and the tem- 
perature in units of the well depth. 

SCOZA provides a closure relation for the OZ equation 
(|l|) by expressing the direct correlation function c(r) in 
terms of the potential v{r), as for other approximations 
like the mean spherical approximation (MSA), PYA, hy- 
pernetted chain (HNC) approximation etc., but intro- 
ducing one or more state dependent parameters that can 
be adjusted to force the system to satisfy various exact 
thermodynamic relations of the system. In particular, 
the simplest SCOZA studied so far assumes 



while c(r) for r > a is composed of two contributions, 
describing respectively the soft part of the potential and 
the hard core, as 

g-6ffs(l — cr) 

c(r) = -APvir) + Khs r > a. (6) 

r 

The Yukawa function in Eq. (^) takes into account the 
contribution to c(r) arising from the hard-core repulsion. 
Thus, the two parameters Kj^s ^nd bns can be deter- 
mined by setting v{r) = in Eq. (^) and requiring that 
both the compressibility and the virial route to thermo- 
dynamics lead to the Carnahan-Starling equation of state 
for a hard-sphere fiuid This amounts to describing 
the hard-sphere correlations via the Waisman parameter- 
ization ||2^. The soft part of the contribution in (j|) is 
assumed to be proportional to v{r) and hence it has the 
same range of the potential. The proportionality con- 
stant A is calculated by imposing the condition that the 
compressibility and the energy routes yield the same re- 
sult. This corresponds to the condition p2[. 
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(7) 



g(r) = 



r < a 



(5) 



where is the excess internal energy . Eq. ^ im- 
plies a partial differential equation (PDE) for (S) . 
The Yukawa potential (^) lends itself particularly well to 
the SCOZA scheme, because for this kind of interaction 
it is possible by means of Eqs. (|^), (^) to establish an 
analytic relation between c(g = 0) and C/*^^, which al- 
lows us to obtain straightforwardly a closed PDE from 
Eq.(^ by using V^'^ as the unknown quantity instead of 
/?) pT| , p^ . Once the internal energy has been ob- 
tained by numerical solution of this PDE, the Helmholtz 
free energy is calculated by integration with respect to 

The reasons which led us to adopt SCOZA can be sum- 
marized as follows. The method uses a reasonable choice 
of functional relationship between c(r) and v(r) on the 
basis of numerous calculations obtained over a number 
of years. Furthermore, it has been shown that the 
results obtained are in excellent agreement with com- 
puter simulations, at least for not too narrow ranges of 
the potential. In this sense it may be viewed as the best 
semi- analytical method to study the Yukawa potential. 
Finally, it is simple and convenient to obtain numerical 
solutions of Eq. , and this is helpful to make a survey 
of a problem in a large parameter window, rather than 
just in a small part of the parameter space. To give the 
reader some sense of the accuracy of SCOZA in some of 
the regimes we will discuss later, we present in Fig. ^and 
^ the phase diagrams for h = 1.8 and 6 = 6 calculated 
using SCOZA, compared with Gibbs Ensemble Monte 
Carlo (GEMC) studies (l6|j2|]. The data by Shukla are 
the largest system sizes with this potential yet studied 
by GEMC. In all cases, the agreement with SCOZA is 
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excellent, and we may, in this regime of 6- values and in 
calculating the phase diagrams, consider SCOZA to be 
equal to the best simulations. The method is superior to 
the other simple closures, including for example MSA. 

Having said this, we note that there is little real infor- 
mation available about the detailed reliability of the clo- 
sure relation where the potential becomes much narrower 
than for 6 = 9, though it is reasonable to suppose that 
many properties are still satisfactory for somewhat larger 
^-values. Another comment we may make is that, when 
the range of the potential narrows, Eq. may not be 
the optimal closure to ensure that c{r) is accurate. Both 
these points should be regarded as words of caution, and 
as potential directions to develop the SCOZA method to 
better represent the structure in such problems. Indeed 
this will be the subject of future work. 



B. The Percus-Yevick closure for the SW potential 

We shall now briefly discuss the closure for a fluid of 
colloidal particles interacting via a square well potential. 



oo r < (7 
vir) = I — e a < r < (7 ^ 
(J + 5 <r 



(8) 



where, in the present discussion, we set e = ct = 1 and 
we define the square well parameter A = (5/((T + 5), 
which parameterizes the attractive range of the poten- 
tial. This model has been already the object of great 
interest in colloidal science 10|. The state of the sys- 
tem is specified by three control parameters, the packing 
fraction = Trpcr^/G ( where p is the number density, i.e. 
p — N/V), the temperature fc^T, and the square- well 
parameter A of the attractive shell. 

The PYA for c(r) is g{r) — Q for r < a and 



c(r) g{r) 



1 



n/9"('-) 



(9) 



outside the hard core l^^l- We solve the OZ equation in 
PYA using Baxter's method of the Wiencr-Hopf factor- 
ization Ip?! , ^ . This correspond to rewriting Eq. (|l|) in 
terms of the real factor function Q(r), defined for r > 0. 
For < r < R, R being the range of the potential {a + 5 
in the present case), one has 



rc{r) = -Q'(r) + 27rp / ds Q'{s)Q{s - r) (10) 
as well as, for r > 0, 



rh{r) = -Q'{r) + 2-Kp \ ds(r- s)h{\r - s\)Q{s) 
Jo 



(11) 



S;' - Q{q)Q{qr 
Q{q) - 



2ttp I dre">''Q{r) 
lo 



(12) 
(13) 



The resulting equations, obtained implementing the PYA 
in Eqs.(|l^) and (^ij), are then solved numerically to cal- 
culate the structure factor Sq (see reference [|l^ for more 
details). 



III. THE PERTURBATION THEORY APPLIED 
TO SOLID AND FLUID PHASES 

A. The Crystal Phase 

In this section we shall discuss the method used to cal- 
culate the solid free energy. The perturbative approach 
has previously been used by Cast et al. [|l5| to construct 
the phase diagram of a colloidal solution with depletion 
interactions. 

Further applications of the method appear in refer- 
ences [p^|Jl6|. For example, in |14| the authors compared 



the results obtained for both phase diagrams and free 
energies of the solid and liquid phases for an Asakura- 
Oosawa potential and made comparisons to Monte Carlo 
calculations. They noted that the phase diagrams are 
overall reasonably satisfactory, but that the crystal free 
energies are excellent, in most regimes being essentially 
quantitative. It was concluded that the potential limi- 
tations in accuracy of the phase diagrams arise from the 
use of perturbation theory to the fluid and liquid phases, 
rather than the crystal. Also, for the crystal free en- 
ergies there does not appear to be a significant loss of 
accuracy when the range of the potential becomes nar- 
row. We have, therefore, applied perturbation theory to 
the crystal state of the Yukawa potential. Besides what 
we show here, it is opportune to note that there are in 
the literature other perturbative approaches, such as in 

ii- 

The method is summarized as follows. We separate out 
explicitly the interaction potential as a hard-core contri- 
bution plus the attractive tail. The hard-core part of the 
potential is used as a reference for the perturbation, and 
the attractive tail is the perturbation itself . In other 
words, we decompose the potential as. 



v{r) = VQ{r) + Vatt{r) 



(14) 



Q(r) determines Sq via its Fourier transform. 



where VQ{r) is the hard-core repulsive potential and ex- 
pand around the reference state VQ{r). With this choice, 
the zeroth order term of the free-energy expansion coin- 
cides with the hard-sphere free energy. Once the pertur- 
bation expansion is carried out to second order we have 
for the Helmholtz free energy the following expression 
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PF PFo Bp f , , , , , ^ 
^ = — + Y y ^att{r)go{r)dr + — 



(15) 



here PF2/N is the second-order perturbation term. Fq 
and go are respectively the Helmholtz free energy and the 
radial distribution for the reference hard-sphere system. 

We focus our attention on the second order term in 
the expansion in Eq. (JTsI). Indeed, its exact evaluation 
requires the calculation of higher order distribution func- 
tions 1^ , which are very hard to compute or to approxi- 
mate reasonably. Barker and Henderson proposed an 
approximation to F2 , based on the following observation. 
Since, 



1 



N 



(16) 



where Wn = J2i<j ^attil^i — Barker and Hender- 
son proposed to divide the space into concentric spheri- 
cal shells and to calculate averaged properties using the 
number of particles in each shell. Following this route, 
they rewrote Eq.^ in terms of average numbers in the 
shells, 



(17) 



where Ni is the number of particles in the shell i and Vi 
is the perturbation energy, considered constant, within 
the shell. The first approximation consists of ignoring 
the correlations between shells, i.e. 



{N,N,)-{N,){N,)^0 



(18) 



for i j. Moreover, inside a given shell, a second ap- 
proximation is made. 



(iVf ) - {N,f « {N,)kBTidp/dP) 



(19) 



The two approximations ( |18| ) and (y^ are equivalent to 
considering the volume of the shells to have the compress- 
ibility properties of a macroscopic portion of space (for 
more details see |32 ). As a result, Barker and Henderson 
approximated the second order term in the expansion as. 



I3F2 
N 



f3p 
4 



dp 
dP 



^=-^1^ vUr)go{r)dr 



(20) 



This approximation was found to be satisfactory in all 
calculations carried out so far. In our work the integrals 
in Eqs. (^5|) and ( pO| ) have been performed by a five- 
point integration rule, while for differentiation a central- 
difference scheme has been used. p^ . 

To carry out the calculation, we require the Helmholtz 
free energy and radial distribution function of the un- 
perturbed hard-sphere system in the solid phase. It has 
been shown by computer simulation that a hard-sphere 
fluid shows a solid-fluid transition, for which the fluid 



phase alone exists up to a packing fraction (j) = 0.49 
and the solid FCC phases exists for cj) > 0.55 |Q. In 
between, there is a two phase coexistence of solid and 
fluid. These properties are well studied and the vari- 
ous information required in perturbation theory can be 
deduced from these studies. We note in passing that 
recently a renewed interest has been arisen in the equi- 
librium structure of a hard-sphere crystal. Indeed it has 
been believed for a long time that hard spheres crystal- 
lize with an FCC structure. Confocal microscopy obser- 
vations, however, have rather found a random hexagonal 
phase which consists of a stack of FCC and hexagonal 
close packing (HCP) layers Simulation seems to ex- 
plain the phenomenon in terms of the small free-energy 
difference between FCC and HCP structures js^] . In this 
paper we assume the crystal equilibrium structure to be 
FCC since this is believed to be the more stable 

To provide continuity with previous authors we make 
the choices described below. The state equation for a 
hard sphere FCC solid has been proposed by Hall |3|] 
who derived a phenomenological expression based on 
computer simulation results, i.e. 



PhsV ^ 
NksT 

- 02 _ o.67825(/)3 - 04 - 0.50^ - 6.0280^/(0) 
1 - 30 302 - 1.00430503 



(21) 
(22) 



with /(0) = exp((7r\/2/6 - 0)[7.9 - 3.9(7rV2/6 - 0)]). 
The compressibility can be derived by differentiating the 
compressibility factor Zhs a-s, 



(23) 



In order to calculate the excess hard-sphere Helmholtz 
free energy Fq^ from the compressibility factor Zhs^ a 
thermodynamic integration in the packing fraction can 
be performed, obtaining 



(5F^- .^^ fiFS' 



N 



N 



) + / (Zhs - 1) 



dr]' 
V 



(24) 



Since the zero-density limit of a FCC crystal cannot be 
represented as easily as the one for a gas, alternative 
routes to perform the thermodynamic integration in ( |2^ ) 
have to be devised We have chosen to perform 

the integration starting from a packing fraction value 
of 0* — 0.544993, for which the value for the free en- 
ergy has been calculated by computer simulation to be 
Fo^"^(0*) N = 5.91889 [||. We recaU that the excess 
Helmholtz free energy is defined as the excess with re- 
spect to the ideal gas contribution [ p7[ . 

For the radial distribution function go{r) we have used 
an analytic formulation proposed by Kincaid and Weis 
that fits Monte Carlo simulation for a hard sphere FCC 



solid |40|. This formulation is known to provide a good 
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estimate of the hard sphere radial distribution function, 
at least in the range 0.52 < cj) < 0.56518. Eqs. ( p^ and 
( pO| ) can now be solved. 

Once the Helmholtz free energy is evaluated, following 
the route we have just described, the Gibbs free energy 
and the pressure can be calculated as 



/3P: 



dp 

pPG pI3F 



(25) 
(26) 



We have earlier noted that the perturbation theory for 
the crystal is highly accurate. Indeed, the second or- 
der perturbation term is useful, but it is interesting to 
note that the great bulk of the free energy correction for 
narrow well problems is captured by the first order term 
alone. To understand this, it is worth reflecting on the 
fact that the free energy of the crystal in the presence of 
short-ranged attractions is referred to the hard core crys- 
tal, and there is no question of the perturbation theory 
having to determine a priori any gross structural infor- 
mation. The corrections from attractions arise by virtue 
of the small changes in local vibrations that the particles 
make around their lattice positions, a portion of these 
motions involving the particles being within their mu- 
tual attractive range. In first order perturbation theory 
these contributions are treated as if the nature and distri- 
bution of the vibrations is unchanged, and the additional 
attractive energy contributions calculated essentially as 
an integral over the attractive potential multiplied by 
the zeroth-order hard core correlation function. The fact 
that second order contributions are typically small is sug- 
gestive. In a heuristic manner, we may argue that the 
intrinsic limitation on the extent and complexity of the 
phase space of the localized particles, imposed by their 
being in a crystalline state, means that even when at- 
traction is incorporated, the changes in the nature and 
distributions of these vibrations is small. We note that 
this rationale is clearly inapplicable to the case of liquid, 
gas, and fluid states, where the addition of attractions to 
significantly affect the distribution of particles motions. 
Possibly this is the reason why perturbation theory works 
well for the crystal, indeed far outside its expected lim- 
itations, but is less successful for the other states. For 
completeness we note that for the crystal (Yukawa and 
SW) the errors, as estimated by the ratio of second-order 
to first-order terms are typically of order 0.5%. For the 
liquid they can be larger. We have however, also studied 
a square-well fluid using perturbation theory, and used 
these results along with those for the crystal to generate 
a phase diagram. The perturbation theory of the square- 
well fluid is therefore briefly discussed below. 



B. Liquid phase for the SW fluid 

In this section we shall discuss the method we adopted 
to calculate the thermodynamical properties of a fluid of 
colloidal particles interacting via a SW potential. 

We chose the hard-sphere fluid as the reference sys- 
tem and treated the attractive part as the perturbation. 
The natural choice to describe the thermodynamics of a 
hard-sphere fluid is the Carnahan-Starling (CS) equation 
of state [Esl , 



I3P l + ri + rf ~rf 



P 



(1 - r^f 



(27) 



The CS equation provides an accurate account of the 
thermodynamic behaviour of the hard-sphere fluid for the 
entire region of the fluid phase. Its very simple analyti- 
cal form makes it possible to obtain a closed expression 
for the Helmholtz free energy by integrating over density, 
as in Eq. (|2^). The zero-density limit of the free energy 
is the ideal gas value, so the thermodynamic integration 
starts from zero density. Thus, we obtain 



N 



^(4-377) 

(1-7^)2 



(28) 



The compressibility is evaluated as in Eq. (^) by dif- 
ferentiation. For the radial distribution function we use 
a modiflcation of the analytical PYA g{){r) for a hard- 
sphere fluid that was proposed by Verlet and Weis ]4l| ] 
to overcome certain limitations of this closure. Indeed, in 
PYA the contact value g^lr — a) of the radial distri- 
bution function underestimates the real value obtained 
by computer simulation and, also, the oscillations of the 
tail are slightly out of phase and too weakly damped. 
Verlet and Weis proposed. 



5o(r/CT, 0) = .go(r-/fT', </-') + ^31 (0 



(29) 



Here, g(r) is written as the sum of two terms. The flrst 
term corresponds to the solution for gQ{r) within PYA, 
but evaluated at a smaller packing fraction value <j}' and 
a smaller diameter cr', while the second term is a short- 
ranged correction 5gi (r) . The parameter a' is then eval- 
uated via a minimization of the difference between the 
simulation result for go{r) and the PYA analytical re- 
sult, between 1.6(7 and Scr. This contribution improves 
the long-range behaviour of the PYA result. The addi- 
tion of the short-range term 5gi (r) improves the value at 
r — a. Analytical forms for (j)' and 5gi{r) arc given by the 
authors pi]] . The improved radial distribution is within 
1% of the computer simulation result in the whole range 
of packing fractions. With this result, the Helmholtz free 
energy for the hard-sphere reference system is calculated, 
and we may then proceed as in the previous section to 
calculate the Gibbs free energy G and pressure P using 
Eqs. (H) and ^ 
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IV. THE CONSTRUCTION OF PHASE 
DIAGRAM OF THE YUKAWA MODEL BY 

HYBRID METHOD 

In the case of the Yukawa potential we have modified 
the approach of previous researchers somewhat in order 
to obtain the benefits of the best methods of condensed 
and hquid state theory. We have used SCOZA to cal- 
culate the liquid, gas and fluid phases free energy, but 
applied perturbation theory for the crystal free energy. 
To mark the difference with previous calculations, where 
phase equilibrium lines had been calculated by perturba- 
tion theory both for the crystal and for the fluid phase, 
we name our approach a hybrid method. 

Phase boundaries between two phases (gas-liquid, 
fluid-solid etc.) are obtained by imposing the standard 
conditions, 

= M^'^ (30) 

p(l) = p(2) (31) 

where n is the chemical potential, i.e. the Gibbs free 
energy per particle, fi = G/N. 

In those cases where SCOZA and perturbation theory 
are quantitatively validated the equilibrium phase dia- 
gram is highly accurate. As noted above, and shown in 
Figures 0and|, the SCOZA is well validated up to values 
of 5 = 9. For values of b less than or equal to 6, results are 
indistinguishable in terms of phase equilibria from the 
best simulations that have been carried out ||2^ , [l6| , ^ . 
Similarly, as we shall discuss later, the perturbation the- 
ory rarely produces an error of more than 0.5% in the 
free energy of the solid phase, although this analysis is 
based on certain assumptions about the perturbation se- 
ries. Combining these observations, we believe that our 
phase diagrams are quantitatively accurate up to at least 
6 = 9. Beyond that, we make no particular claim, except 
that we expect that this hybrid method should still re- 
main superior to the typical theoretical approximations 
that have been applied previously. Simulations have not 
been carried out beyond 6 = 9. 

V. GLASS TRANSITION AND MODE 

COUPLING THEORY 

A. Theory 

The study of the glass transition in colloidal systems 
has been one of the most striking cases of veriflcation 
of the current theories of super-cooled liquids. Early 
experimental studies |^,^ involved colloidal particles 
that are very closely represented by hard spheres where 
only excluded volume effects are important at high con- 
centrations. Moreover, in contrast to simple atomic liq- 
uids, it is possible to avoid the crystalline phase beyond 



volume fractions of 49% for sufficiently long periods of 
time to study the glassy-type dynamical processes, and 
ultimately the colloidal glass. The agreement between 
certain aspects of MCT and experiments on colloids is 
quite satisfactory and the details of the time corre- 
lation functions are quite well reproduced. It is widely 
believed that, in deeply supercooled molecular liquids, 
the slow dynamics involves more complex dynamical pro- 
cesses than those described by MCT, and there the the- 
ory becomes of more qualitative applicability. Thus, the 
case of colloidal particles is of some practical interest in 
applying this type of theory. 

In fact, even for colloids, small discrepancies appear 
in the comparison between experiments and MCT. The 
most important is the value of the critical volume fraction 
for the hard-sphere arrest transition, the experimental 
value being about 58%, while MCT predicts about 52%. 
This is of little importance where the dynamical laws at 
the hard-sphere transition are being compared between 
experiment and theory. Previous researchers have ap- 
plied a shift to the transition volume fraction, and then 
fitted the laws in this region . Since the only current 
information on arrest driven by attractive interactions is 
that provided by MCT there is as yet no accepted 
manner in which we can correct the MCT curves. This 
is somewhat inconvenient in the current context since for 
some parts of the parameter space the equilibrium phase- 
diagrams are quantitatively accurate, and it would be 
very satisfying to be able to superimpose, without cor- 
rection, the relevant MCT arrest curves. 

We now briefly review the nature of MCT, and discuss 
the type of information it yields. The MCT of super- 
cooled liquids describes the non-ergodicity transition by a 
nonlinear integro-differential system of equations for the 
normalized time correlation functions of density fluctua- 
tions ^(q,t). Apart from parameters entering from the 
microscopic motion, the only input to the MCT equa- 
tions is the equilibrium wave-vector dependent structure 
factor of the system, Sq. The glass transition lines can 
be identified by studying the long time limit of the MCT 
equations, which determine the non-ergodicity parame- 
ter of the system fq = limt^oo ^{q,t)- An ergodic state 
is characterised by /q = 0. This value is always a solu- 
tion of the MCT long-time limit equations Thus, 
the glass transition appears as an ergodic to non-ergodic 
transition for the system, where fq ^ solutions arise. 
These points, thus, correspond to bifurcation singulari- 
ties of the MCT equations, and, depending on the num- 
ber of control parameters of the model, these can be of 
increasingly higher order, producing interesting features 
of the arrested states diagrams. 

A good quality Sq is an important input for a good 
description of the MCT arrest transition, as for the equi- 
librium phase diagram. In the earliest discussions of col- 
loidal systems with short-ranged attractive interactions 
the Baxter interaction [B3,E6|, a limiting case with an 
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infinitely deep and zero-ranged SW attractive potential, 
was discussed by a number of authors Sub- 
sequent studies indicated that the MCT equations are 
pathological for this interaction [|9|. The calculations 
were therefore extended to a SW potential both in the 
PYA and the MSA |o|. Another solution of the MCT 
equations was obtained using the Yukawa potential and 
the MSA jl9|. Some common aspects emerged in these 
works. However, the SW model at first appeared to give 
a richer behaviour for the arrest transition curve. Thus, 
the system was shown to possess a glass transition curve 
in the parameter plane (0, ksT) the shape of this curve 
depending on the value of the SW parameter A |]lo| . For 
narrow well-widths, two branches of the glass curve have 
been identified. These have been interpreted respectively 
as a transition between a fluid phase and repulsion dom- 
inated glass (this is the typical repulsive glass) and a 
transition between a liquid and an attractive-interaction 
dominated glass (named the attractive glass). The two 
branches join and for approximately A < 4.11%, a glass- 
glass coexistence between the two different types of glass 
appears. This coexistence line terminates in an end- 
point, beyond which the non-ergodicity parameters be- 
come the same for the two types of structures. The rele- 
vant singularity points, such as the end-point or the point 
where the glass-glass transition line reduces to a single 
point, are identified with higher order singularities of 
MCT equations and lead to unusual logarithmic dynami- 
cal relaxation laws |l^ jl^ . The mechanical properties of 
the system have been also studied |l^ and they reinforce 
this interpretation. Earlier studies of the Yukawa poten- 
tial did not locate this glass-glass phenomenon |jl^ , but it 
was subsequently realized that the screening parameters 
which had been studied were not large enough | ]50[ |, and 
further calculations seem to give clear indications that 
both the SW and Yukawa potentials give the same typi- 
cal behaviour ||5l| , implying that this does not crucially 
depend either on the potential shape or the approxima- 
tion used for calculating the structure factor. Thus, it is 
now believed that this glass-glass scenario, and the at- 
tendant dynamical laws, is essentially a universal feature 
of the very short-ranged attractive potential. 

We note in passing that the formation of two solid glass 
phases for very short-ranged potentials should not be too 
surprising. In fact we have earlier alluded to the fact that 
there are two crystalline phases in the phase diagram of 
such potentials. We may typically view glasses or ar- 
rested states as long-lived meta-stable states of the sys- 
tem that have not been able to equilibrate to the nearby 
crystal, and that are trapped in a restricted portion of 
phase space. In this sense we may expect each crystal to 
have associated to it a particular glass-type. Since one 
of the crystals in our phase diagram is "attractive en- 
ergy dominated" , and the other "repulsive energy dom- 
inated" , it is hardly surprising that there should be two 
types of glass, dominated by the two regimes of interac- 



tion. 

A further comment on the relevance of these glass 
curves is that we may see them as more than simply 'tran- 
sition curves'. If one reflects more deeply on the nature of 
the equilibrium phase diagram, and the thermodynamic 
states present in them, we recognize that at a deeper 
level they are reflecting the fact that for those particular 
parameters, the phase-space is dominated by a particu- 
lar structure: a crystalline, liquid, or gas structure. The 
arrest curves carry analogous information. Thus, in the 
vicinity of the arrest curve, we may understand that most 
of phase space is becoming increasingly inaccessible, and 
breaking into smaller regions that are disconnected. That 
this may occur even for states that appear to have static 
structures typical of liquids is the distinguishing feature 
of glasses. Associated to this observation are those of dy- 
namical slowing and tendency to arrest during any phase- 
separation through which such a glass curve passes. We 
shall argue in the conclusions that such phenomena are 
relevant in protein crystallization. 

B. Experimental Studies 

We now discuss some of the particular experimentally 
determined features that are associated with dynamical 
arrest driven by attractive interactions. For example, 
Verduin and Dhont Q determined a curve of structurally 
arrested states in the phase-diagram of a system with 
short-ranged attractive depletion interactions. This lo- 
cus, in some cases, intersects the binodal line and is re- 
ferred to as transient gelation when observed in the spin- 
odal region. Interestingly, it was these experimental Q| 
authors who were the first to comment that MCT might 
also be applicable to cases where attractive interactions 
are important. As we noted above, subsequent results of 
such calculations have been most interesting ||lo| , [ll| , |l9| . 
In addition a number of other experimental programs 
involving particles with depletion interactions have been 
published which offer many interesting insights for exam- 
ple in ||5^ , |53|] . In particular, Poon et al [|| have studied 
the arrest transition for systems where the range of the 
interactions is short, and their more recent work in this 
topic involves detailed connection to the theory described 
above |5j]. Other systems may have certain advantages 
over the depletion interaction system, but it is as yet too 
early to decide this issue [|j5|-|^. The results are typ- 
ically quite promising, with some of these other systems 
also exhibiting some of the phenomena predicted by the 
theory. We may note in particular a more recent set of ex- 
periments that are intriguing in that they make detailed 
predictions for the correlation functions in a particular 
(re-entrant) part of the phase diagram. Thus, concentra- 
tion time correlation functions have been observed in a 
polymer micelle system with a decay process much longer 
than the usual stretched exponential, and the results have 
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been well- fitted to a logarithmic time- relaxation, |61|, as 
predicted by the theory |^ . 

The development of experimental understanding, and 
deepening of the theory of systems with short ranged 
potentials is really just beginning, and many experimen- 
tal programs have now been commenced or reoriented 
in efforts to make progress. However early information 
indicates that the MCT type theory may be able to de- 
scribe main elements of the principal phenomena, at least 
for reasonably high volume fractions, where it is possible 
to separate aging from dynamical arrest in a reasonably 
clean manner. 



VI. ARE SPINODALS IN COLLOIDAL SYSTEMS 
MEANINGFUL? 

Here we will take the liberty of raising a few issues 
in relation to spinodal curves that are calculated via 
SCOZA, or indeed many other typical liquid state the- 
ories. The reason that we make these comments is that 
such curves should have a particular status for these col- 
loidal systems that are not relevant generally for molec- 
ular systems. We note first that the spinodal curve is 
determined from the condition that the curvature of the 
free-energy with respect to the relevant density variable 
becomes zero, and that this corresponds to the fluid 
phase becoming unstable as we lower the temperature. 
Between the binodal and the spinodal, the liquid and gas 
states in co-existence may be the global free-energy min- 
ima, but the fluid state remains meta-stable. Inside the 
region bounded by the spinodal, only the liquid and gas 
states in coexistence are stable. Now it is well known |6^] 
that the free energy is a convex function, and it possesses 
only one minimum, and for some years now it has been 
understood that the spinodal curve determined from ap- 
proximate theories (eg. mean-field theories) that con- 
sider two separate branches of the free energy and then 
connect them, has no real scientific basis. Careful Monte 
Carlo simulations carried out in systems of increasing size 
have lead to the conclusion that the spinodal curve 
shifts with the system size, merging with the binodal in 
the limit of infinite systems. There is some loosely de- 
fined kinetic phenomenon however |l6^ , |65| , though even 
there it is not possible to define a spinodal curve, but 
a cross-over regime where the kinetic mechanisms begin 
to change from nucleation and growth to more collective 
phenomena. Interestingly enough, when the particle size 
becomes large these more sophisticated expectations are 
less relevant. Thus, it transpires that the relevant pa- 
rameter in this story is the ratio of the particle diameter 
to the correlation length of the fluid. For very large par- 
ticles such as high molecular weight polymers, colloidal 
and other particles, the microscopic length is so large 
that one has to be extremely close to the critical point 
to see fully developed fluctuations beyond the mean-field 



type ideas. Another consequence of this is that critical 
exponents in such systems as proteins [ |66| and micelles 
p7t have often been measured with mean-field values be- 
cause experiments were not performed in the true critical 
regime. Similarly, the normal scepticism about the exis- 
tence of a spinodal curve should be less relevant here, and 
we may expect the colloidal systems to exhibit quite rea- 
sonable spinodal behaviour. We have therefore included 
the spinodal curve in our phase diagrams. 
In concluding this section on dynamically defined objects 
in the phase diagrams, it is commented that we have 
chosen to plot the MCT curves through the meta-stable 
regions between binodals and spinodals. Again, in col- 
loidal systems, for the reasons given above, it is to be 
expected that such curves would have meaning, whereas 
they would not be meaningful in molecular systems. 



VII. PHASE-EQUILIBRIA AND DYNAMICAL 
ARREST LINES 

In this section we combine the results from the differ- 
ent techniques described above to exhibit the state of the 
system for given well depth. 



A. Yukawa potential 

In this section results for the Yukawa potential with 
the temperature in units of the well depth (fcsT/e) plot- 
ted against the volume fraction of the system, are pre- 
sented. The values of the well width are determined via 
the screening parameter of the Yukawa b. The hard core 
radius is fixed at unity, so all quoted lengths are in units 
of the hard core diameter. We show results for values of 
6 = 5 (Fig. I), 6 = 6.05 (Fig. |), 6 = 7.5 (Fig. 5) 6 = 30 
(Fig. I), 6 = 60 (Fig. I), and b = 100 (Fig. |). For 
comparison we also show (Fig. ^) a calculation for the 
square well system with A = 0.03. 

We begin with the largest well width, corresponding to 
6 = 5, see Fig. ^. Here the well width, considered say as 
the distance of half amplitude of a Yukawa, is comparable 
to the particle hard-core size. This is the typical situation 
that we are familiar with in elementary phase diagrams 
of atoms and molecules where van der Waals interactions 
predominate. Thus, we see the expected pattern of phase 
behaviour. Below the critical temperature, the gas-liquid 
phase equilibrium occupies the greater part of the low- 
and middle-range of densities, above the triple point. The 
crystal is favoured at higher density, and the liquid- and 
fiuid-crystal boundary is nearly vertical, that is at fixed 
density, refiecting the substantial absence of any energy 
scale in the problem. The crystal is so tightly packed, 
and the attractions are so spread out across the system 
that it is only the repulsive part of the potential that is 
fundamental for crystallization. Indeed this is one of the 
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important ideas in traditional liquid state theory; that 
attractions are not relevant to crystallization. 

The asymptotic limits of solidification and melting 
boundaries of this coexistence at high temperatures re- 
flect the hard-sphere limits respectively of </) ~ 0.49 and 
(j) ~ 0.55, as expected from many simulations and theo- 
retical observations [ p"5||34| ] . The triple point temperature 
is labelled Tp. Within the crystal phase there is also the 
boundary for the dynamically arrested state (the MCT 
transition line), again with high-temperature asymptote 
of (/) ~ 0.52(circles), the hard-sphere volume fraction for 
the glass transition predicted by MCT. We recall that 
MCT underestimates the glass transition packing frac- 
tion by about 0.06. To call attention on this shift we 
report the true MCT curve (circles) as well as the MCT 
curve shifted by 0.06 in packing fraction using * as sym- 
bol. This boundary is also almost vertical, again re- 
flecting the fact that, for wider wells, the arrest tran- 
sition is driven essentially, at high enough temperature, 
by the repulsive part of the potential. Thus, no attrac- 
tive glass is observed for this range of the potential. It 
has been shown theoretically , experimentally ||6|] 

and by simulation that on decreasing the range of 
the attractive potential the fluid-fluid coexistence curve 
becomes meta-stable with respect to the fluid-solid one. 
This means that for short enough potential ranges the 
liquid-gas critical point is hidden in a phase separat- 
ing region. This has two main consequences. The co- 
existence curve between the fluid and the solid is very 
broad implying, for low temperatures, coexistence be- 
tween a very low density fluid and a high density solid 
(FCC). Also, the critical point becomes meta-stable with 
respect to gas-solid co-existence. We will return to some 
of the possible implications of all this later in the dis- 
cussion. We have calculated the particular value of the 
screening parameter, 6, at which the critical point be- 
comes meta-stable with respect to the low density fluid- 
gas equilibrium, to be b* — 6.05, and then presented the 
phase diagram at that value (see Fig. ^). We note that 
slightly different values have been previously reported in 
the literature. Thus, Hagen and Frenkel quote the val- 
ues b* — 7.4 using a Monte Carlo perturbation theory, 
and the value b* — 6 based on Gibbs Ensemble Monte 
Carlo (GEMC) Their GEMC value is very close to 
the value we have found. However a number of other 
values have also been quoted in the literature. Menderos 
and Navascues []69| used a density functional approach 
to determine b* = 8.25, whilst Shukla quotes a much 
higher value of &* — 13. This latter value in particular is 
much higher than previous ones and our value, and con- 
sidering that it is based on more extensive simulations 
than previous research this might be a matter of con- 
cern. The issues in relation to the accuracy of SCOZA 
have been addressed for modest values of b in Figures 1 
and 2, where the liquid-gas phase diagrams were shown 
for 6 = 1.8 and 6 = 6. However, if we examine Table 4 



and Figure 4 of [g8| , we can readily compare the simula- 
tions to the precise predictions from SCOZA for selected 
values of the vapour and liquid densities for a range of 
temperatures. For b* = 6, we find remarkable agree- 
ment between SCOZA and the Shukla gas-liquid equilib- 
ria. Also, whilst they are not quoted, we may estimate 
the critical temperature for differing screening parame- 
ters from these same simulations, and interpolation of 
others, and again conclude that the discrepancy between 
what we find for 6*, and the Shukla result cannot arise 
because of a difference between the SCOZA and GEMC 
results for gas-liquid systems. The problem lies in the 
estimation of the crystal free-energy, or entropy, or esti- 
mation of the triple point, the other aspect of what we 
need to know to determine b* . In the freezing tran- 
sition location is determined by the so-called one-phase 
entropic condition |]70| . This condition implies that the 
freezing density is essentially constant, and almost un- 
changed for the range of b values 1.8 to 10 In fact, 
in Figure ^, the phase diagram for b* = 6.05, we do see 
rather significant deviation of the freezing density in the 
vicinity of the critical temperature from its high temper- 
ature limit. Indeed, comparing the other phase diagrams. 
Figures 5-^, we conclude that this variation is intrinsic to 
the whole short-ranged scenario, since it is the prelude 
to the splitting of the solid phase into two crystalline 
phases. The true underlying discrepancy arises because 
the entropy criterion implies that attractive forces are 
not important in the regime of crystallization currently 
under discussion, whereas the perturbation theory im- 
plies that they are highly significant. Both approaches 
are approximations, but it is possible to estimate the er- 
rors arising from the perturbation theory by considering 
the higher order terms. 

Thus, for b* = 6.05, (here kBT*/e = 0.454 and 
(/>* = 0.230) we can propose to bound the errors in the 
perturbation theory of the crystal by reporting the ratios 
of the second- to first-order terms in the perturbation se- 
ries. In the regime where the low-density fluid and crystal 
are in equilibrium we find that the ratio of the second- 
to first-order terms is approximately 0.0044, whilst the 
first order-term has an absolute value of 10.4 [0. If this 
ratio represents a true estimate of the errors, then the 
perturbation theory would appear to be quite satisfac- 
tory, the curvature of the freezing curve genuine, and the 
attendant fiattening of the fluid side of the coexistence 
also quite accurate. 

In general by estimating the impact of such errors in 
perturbation theory, we can estimate the shift of the fluid 
side of coexistence, and thereby estimate errors. In so do- 
ing we find the typical error in 6* will be less than 1% 
from this source. If we accept this means of character- 
izing the error in perturbation theory, a matter clearly 
based on the series being well-behaved at successive or- 
ders, then we might conclude that the value of b* = 6.05, 
in agreement with Frenkel's original calculation, is a good 
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approximation. If this were confirmed, then the calcu- 
lations reported here for the overall phase diagram are 
probably amongst the most accurate for modest values 
of the range parameter, despite the fact that simulation 
is never used. This is not the primary motivation of our 
paper, but it would be an interesting way of approaching 
phase diagrams in future. Despite these optimistic esti- 
mates, further more careful evaluations are required by 
different methods to find a truly accurate value. 

We should not imply that the value of b* is of such 
crucial importance in the overall picture offered here. 
However, it does provide a useful check between different 
researchers and methods of approximating the phase be- 
haviour in that its accurate estimation requires some sat- 
isfactory and simultaneous treatment of gas, liquid and 
solid phases. 

Finally, we note that ten Wolde and Frenkel have made 
some interesting comments in relation to the kinetic pro- 



cesses that might be expected in this regime |72 



In Fig. 5 the case b — 7.5 has been plotted. As noted 
by a number of authors, the meta-stable gas-liquid phase 
equilibrium curve has now flattened considerably, and the 
low-density fluid-crystal coexistence (on the fluid side) 
has nearly the same slope. Increasing curvature of the 
crystal side of the fluid-crystal equilibrium curve is ob- 
served, arising from the increased influence of attrac- 
tions on the crystal. The binodal line and the spin- 
odal line in this case lie completely within the region 
of fluid-solid phase separation, the triple point has dis- 
appeared and the critical point is buried below the flat 
part of coexistence curve between the low-density fluid 
and the solid. As we mentioned earlier, this meta-stable 
behaviour would normally not be observable due to fluc- 
tuations, but for colloidal systems, globular proteins and 
other large particles, we may expect to observe such phe- 
nomena. Thus, on quenching such a system we might 
expect to see rather a rich pattern of behaviour, depend- 
ing on the density that we quench at, and the depth of the 
quench. In particular, it is noted that we should be able 
to see a meta-stable gas and a liquid, the latter arresting 
into a glass at sufficiently low temperature, because the 
glass curve crosses the binodal, and spinodal at a finite 
temperature. We have earlier alluded to the idea that 
critical fluctuations can play an important role in the for- 
mation of crystals, for example protein crystals f7% . The 
present screening parameter regime would exemplify this 
type of phenomenon since here we have a meta-stable 
liquid and gas that are critical (for large particles this 
should have some observable lifetime). The equilibrium 
phase diagram exhibits a fluid-crystal co-existence, so we 
have the possibility of crystal nucleation and growth phe- 
nomena in the presence of this meta-stable critical fluid, 
and it is this matter that ten Wolde and Frenkel have 
discussed ||7^]. However, we also note the point that for 
this value of the screening parameter, as yet, the glass 
curve has not begun to severely interfere with the gas- 



solid equilibrium curve, and this is an additional advan- 
tage in the formation of crystals. This is in distinction to 
subsequent phase diagrams where the glass curve extends 
across much of the space. 

It is in fact worth reflecting on the shifted, more realis- 
tic placement of the glass transition curve (stars). Thus, 
we see that at high temperature, as expected, from 55% 
to 58% volume fraction values, we have a crystalline state 
that is not interrupted by a glass transition. However, as 
the temperature is lowered, the increased importance of 
the attraction leads to the glass curve crossing the low 
density fluid-crystal coexistence region, and beneath this, 
the crystal may never form without the glass being an al- 
ternative long lived state. This comment is relevant also 
to the case oi b — 6.05, but for b ~ 7.5, the curvature 
of the solidification curve has increased greatly, so this 
effect is emphasized. This dramatic interruption of the 
crystallization scenario will become more and more sig- 
nificant as the range of the potential narrows, and will 
be an important theme in our discussion. 

Now we turn our investigation to the case of very nar- 
row potential ranges. In Fig. || we present the case 
6 = 30. The two-phase coexistence of fluid and solid 
now occupies a much larger portion of the parameter 
space. The gas-liquid critical point is hidden well below 
the fluid-solid coexistence line and on the crystal-fluid 
phase boundary the effects of critical fluctuations will be 
much less. It is interesting to note that the crystal side of 
the coexistence curve exhibits a strong deviation towards 
higher densities. Thus the coexisting solid will be much 
more dense that in the previous cases, since the strong 
short ranged attraction is pulling the particles closer at 
low temperature. This phenomenon is the precursor of 
a solid-solid phase coexistence that in this case is still 
meta-stable and lies within the sublimation curve. 

The glass transition curve is most interesting. Of 
course, for T oo, the curve asymptotes to the hard 
core value. The attractive forces at low temperature now 
begin to strongly affect the curve so that it now turns 
sharply to the left, passing very close to the submerged 
critical point. Thus, we may tentatively assign b** — 30.0 
as that value of the screening parameter at which the sub- 
merged critical point becomes submerged by the glass, as 
well as by the crystal-gas curve. The inherent inaccura- 
cies in the MCT estimate, alluded to before, may mean 
that the b** value may not be very accurate. However 
the phenomenon is interesting. It means that the meta- 
stable critical fluid is now competing with a glass tran- 
sition. It will transpire that the nucleation rate of these 
fluid-crystal equilibria are very low, possibly due to the 
high interfacial tension and therefore it is feasible that 
one may be able to approach the liquid gas equilibrium 
and its nearby glass transition without signiflcant inter- 
ference of the crystal. The type of slowing that would 
arise from a combination of critical slowing down and 
glassy slowing down has not been discussed in the lit- 
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erature previously. It would be an interesting problem. 
However, the phenomenon discussed here would signifi- 
cantly affect the possibility to form high-quality crystals, 
perhaps rendering it essentially impossible. Even though 
the glass may eventually decay in favour of the crystal, 
this will never lead to high quality crystals. Of course, 
an alternative view of this situation is that under these 
conditions it may be possible to make interesting mate- 
rials that have critical fluctuations frozen into the glassy 
phase. Protein scientists wish to make good quality crys- 
tal, materials scientists often wish to make interesting 
materials. Our comments are applicable to both situa- 
tions. 

It is interesting also to note that in this case the 
glass curve passes close to the gas-liquid critical temper- 
ature, but then dips and intersects the binodal line below 
the critical density, passing through the spinodal region. 
Such a scenario has been found by Verduin and Dhont 
Q in experiments on colloidal systems. We should stress 
that for low densities the glass curve itself may not be 
reliable, as was discussed in ||l^. The situation in rela- 
tion to this point is, as yet, not settled. 
We now discuss the case of extremely narrow wells. It is 
worth noting that the previous two values of the screening 
parameters represent the typical range of values accessed 
by those studying depletion-induced attraction between 
colloidal particles Q or globular proteins |l|,^ . The next 
set of phase diagrams (Fig.^ and Fig.|^ represent the 
limit of these types of interaction, and may correspond 
to cases such as grafted coatings on latex particles where 
we can access much narrower ranges of potentials. 
In Fig. 1^ and Fig. || we present the cases 5 = 60 and 
b — 100 respectively. For numerical reasons, it is diffi- 
cult to extend fluid-solid coexistence curves in the low- 
density region of the phase diagram. This is evident from 
the truncation of the phase diagrams, and has no fun- 
damental significance. As already noted in this paper, 
SCOZA has not been tested with simulations for such 
narrow ranges of the attractive potential. On the other 
hand, the main phenomena appear using other methods, 
so at least these are expected to be reliable. We return 
to this point at the end of this section. 
In both figures we observe similar features. The spinodal 
line is now buried deep in the sublimation curve as is the 
critical point. A most striking feature of these phase di- 
agrams is the presence of a solid-solid coexistence. This 
first-order phase transition was already present for less 
narrow ranges but there it was meta-stable (see for exam- 
ple the shoulder in the phase diagram for b — 30). This 
phase boundary represents the coexistence of two crys- 
tals with the same lattice structure but different lattice 
spacing and consequently different density . It is ter- 
minated by a critical point of the solid-solid coexistence. 
The origin of this co-existence is interesting. The pres- 
ence of short-ranged attractive interactions causes com- 
petition with the hard-core repulsive interaction. The 



fact that both have variations that occur on very short 
length scales means that the system may be forced to 
'choose' between the attractive-dominated crystal and 
the repulsive-dominated crystal. 

In the low-density crystal region, i.e. (f> < 0.65, the crys- 
tallization is dominated by entropic effect. In other words 
the the system chooses to optimize the entropy to form 
an FCC structure. Increasing the density the particles 
become closer and at some density they are forced to re- 
main in the attractive shell of their nearest neighbours. 
When this happens, there is a decrease in energy which 
leads to an "attractive" crystal. It is the energy that 
stabilizes the phase. It is indeed clear that such a phe- 
nomenon can be present only if the range of the attrac- 
tive potential is short enough. We may note that an 
isostructural phase transition has been already discussed 
theoretically for other kinds of potential characterized by 
a short-range potential and indeed was also detected 
by simulation 1 12 ij] . It is almost certainly a genuine 
phenomenon. Here the isostructural phase transition is 
present for both 6 = 60 and b = 100. Decreasing the 
range moves the critical point of the transition to higher 
density, and indeed this is also in agreement with the be- 
haviour in simulations [|l2|. It is interesting to note that 
it is possible to find a triple point Tp at which the two 
solids and the fluid coexist at the same temperature. 
In Fig. ^and Fig. jsjthe glass fines have been also plotted. 
They both tend to the hard-sphere limit for high temper- 
atures, and bend towards low densities with decreasing 
temperature, as we have seen in the previous cases. For 
short enough interaction ranges the glass transition line 
does not pass close to the critical point. We note that 
for b = 100 the glass curve appears to break into two 
branches, with an apparent discontinuity at that point 
marked Pd in the figure. The low-density branch is called 
attractive glass, whilst the right-hand branch is called 
repulsive glass. For very short ranges, in other types of 
attractive potentials, we have located a glass-glass tran- 
sition, a transition between two different type of glasses 
originated either by repulsion or by attraction [ |lO[ . For 
the Yukawa fluid such a phenomenon also appears to be 
present, although it has not yet been investigated in de- 
tail |0|j5|. 

Both the results for the Yukawa potential considered here 
and those for the square-well fluid that we present next 
for comparison clearly show that the distinction between 
the attractive solid and repulsive solid becomes sharper 
as the range of the potential becomes narrower. 
To present in a coherent way the role of the attractions 
on the crystal, glass and liquid-gas coexistence line, we 
show in Fig. ^ the dependence of Tc on b and both 
the glass transition temperature and the solid-fluid first- 
order transition temperature at the critical packing frac- 
tion (/)c. For completeness, the inset shows 0c (^)- 
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B. Square- well potential 

Finally we discuss a single example for the SW model. 
By so doing we wish to make the point that the main phe- 
nomena that have been discussed above are independent 
of the details of the shape of the potential, and are essen- 
tially universal. We note, however, that the SW phase 
diagram is not expected to be so quantitatively accurate 
as that of the Yukawa potential for reason discussed in 
section 



III 



The SW model was solved as discussed in section [IB 



We discuss a case where the range of the potential is 
very narrow, i.e. A = 0.03, and the result is presented 
in Fig.^. It is clear that the situation is very similar 
to that for the Yukawa fluid, i.e. a solid-fluid phase co- 
existence extends from high temperatures (where again 
it reaches the correct hard sphere limit) expanding dra- 
matically towards low and high densities for low enough 
temperatures. An isostructural solid-solid phase tran- 
sition, with a critical point, is also present in this case. 
Indeed a limited correspondence between different poten- 
tials based on their general characteristics (effective core, 
range and energy scale) has been recently proposed , 
and it is possible that this idea may have more general 
applicability. We hope to return to this more general 
concept of corresponding states at a later point |[75| . 
The glass line also has a very similar shape to the Yukawa 
fluid. As noted above, and in earlier publications, we 
have a glass-glass transition, terminating in a higher or- 
der glass singularity, the A3 transition |]10[. We believe 
that the presence of these two types of glass, the attrac- 
tive and repulsive glass, is the disordered analogue of the 
presence of the two types of crystal discussed in some 
detail above. In this case we may draw also an analogy 
between the presence of the isostructural critical point 
and the presence of this MCT singularity A3 end-point 
at which the two glasses become identical. Also as the 
well width gets larger, the crystal-crystal critical point 
vanishes. We believe that the glass-glass analogous of 
this is the A4 point [0[l8|. 

It is quite reasonable to suppose that for every crystal, 
there should be an analogous glass, and for every critical 
point of such a crystal-crystal equilibrium there should 
be such an MCT singularity. It will be interesting to 
explore this idea in future. It has the appeal of a poten- 
tial general joint classification of equilibrium and glass 
transitions. 



VIII. CONCLUSIONS 

We shall use the conclusions section of this paper for 
two purposes. Thus we shall attempt to sum up the 
practical conclusions of our calculations, but, at the 
same time, try to considerably broaden the discussion 



to make contact with the main experimental situations 
where they might be useful. So far, we have focused the 
discussion quite strongly on the narrowly defined conse- 
quences of studying a short-ranged hard core attractive 
Yukawa potential, so we will begin by summing up that 
aspect of the discussion. 

Firstly, from the technical point of view, we have achieved 
a certain success combining a good liquid-state method 
of calculation with perturbation theory. The resulting 
hybrid method takes free-energies from SCOZA for the 
gas and liquid states, and from perturbation theory for 
the crystal. We have indicated, in broad terms, how this 
overall strategy could be applied to phase diagrams in 
general, and how it may be qualified by checking of er- 
rors and relating these errors to shifts in the phase bound- 
aries. Given the potential to exploit powerful methods 
of liquid-state theory, and the remarkable success of per- 
turbation theory for the ordered state when one chooses 
the correct zeroth-order state, this may be a competitive 
manner in which to proceed for many problems in future. 
On the other hand, the methods to determine dynamical 
arrest (eg. MCT) are not nearly so developed, despite 
their relative success in colloidal science. In particular, 
the absolute values of the density and temperature at 
which the arrest takes place is not correct. This should 
not be surprising. Equilibrium theory has had the ben- 
efit of many more years of development, and much more 
effort devoted to bring it to this level of achievement 
whereas for the dynamical arrest alternative routes or 
different approaches have been developed only recently 
| |76[ . However, this aspect is quite inconvenient. 
As we try to apply this theory to more realistic situa- 
tions in colloids, materials and biology, we see increas- 
ingly this important motif of a competition between the 
equilibrium Boltzmann view of the matter, and the dy- 
namically arrested aspects. So far, these two fields have 
developed somewhat in independent manners. However 
from the point of view of these practical topics in na- 
ture, there is no distinction, and it is often not the sep- 
arate behaviours, but the interplay and competition be- 
tween them that is primary to the scientific issue. This 
paper is one of the first attempts to make connections 
between these phenomena, but one can clearly see the 
limitations. For greater insight into these topics, we will 
have to address the possibility to both improve system- 
atically the methods of studying arrest transitions, and 
also their consistency with equilibrium transitions. This 
must represent one of the important technical challenges 
in coming years. 

From a broader perspective, we have shown that when 
the range of the potential becomes short in comparison 
to the core size, the subtle interplay between entropy 
and energy begins to change its character. The range of 
densities over which configurational entropy is relevant 
is much reduced, and one begins to lose the liquid state, 
in favour of crystals, or arrested glassy states. The rea- 



14 



sons have been discussed at the beginning of the paper. 
In essence they amount to the fact that to retain the 
benefits of short-ranged attractions the particles must 
not depart too much from their typical inter-particle dis- 
tance, or they are no longer in their mutual attractive 
well. This loss of freedom of motion, and restriction of 
favourable configurations, leads to a lowered configura- 
tional entropy. Another way of expressing the same idea 
is that the short-ranged potential leads to the loss of 
easy fluctuations that can open the cage of neighbouring 
particles that trap a central particle. The probability of 
finding such an 'opening' of the cage is much reduced, 
and the time during which a particle is localized by its 
neighbours increases, diverging at the arrest transition 
to form the 'attractive glass' that we have discussed. As 
the range of the potential narrows, the means of egress 
permitted to the particles is further limited, and the at- 
tractive glass becomes more favoured. This glass is there- 
fore an effective competitor to the liquid and crystalline 
phases of the system, and this is reflected in the fact that 
the fluid phase is eventually erased by the glassy phase, 
and the critical point submerged underneath the curve 
of arrest transitions. This aspect should not be confused 
with the equilibrium phase diagram, although it is inter- 
esting that there are many parallels between the two. 
As we discussed before, when the range of the poten- 
tial becomes very short, the competition between entropy 
and energy is responsible for the formation of two distinct 
crystal phases. Coming from the fluid side, a crystalline 
phase dominated by repulsion is present. This state arises 
from the fact that the, at such densities, the entropy of a 
system made of particles free to move only in their own 
Weigner-Seitz cell of a FCC structure is larger then that 
of the metastable fluid . We named this phase a repul- 
sive crystal by analogy with the glass. On increasing the 
density, this crystal becomes unstable and makes a tran- 
sition to a smaller Weigner-Seitz cell, so that particles 
mainly stay within the attractive well. The structure of 
this phase, that we named attractive crystal, is a more 
compact FCC crystal that optimizes the free energy by 
means of reducing the potential energy. The result at 
this point is that the crystal structure may adopt two 
different states. The situation is analogous to the case of 
a gas-liquid phase separation where the gas (low density) 
phase optimizes the entropy and the liquid (high density) 
has a much lower energy. 

The implications of all of this are profound for practi- 
cal situations. For reasons given earlier, many systems 
such as colloids, globular proteins, "buckyballs" , nano- 
particles, pre-ceramic particulates, and others have this 
property of a short-ranged attraction. In all these cases 
"precipitation" , "gellation" , "glassification" , or solidifi- 
cation are frequently the commonly observed outcome. 
In cases where we consciously seek to make such a state 
this is satisfactory, and it remains only to adjust the po- 
tential to have sufficiently short range to obtain the re- 



quired properties of the solid. However, for cases such as 
globular proteins, and nano-scale or meso-scale ordered 
materials with prescribed optical properties the situation 
is quite different. Here we seek to make a crystal. In fact, 
reviewing the phase diagrams in Figs. 5-|^, we can see 
why the crystal is hard to access. If we work to the right 
side of the glass curve (point marked A in Fig. ||), uni- 
versally we may expect to fall into the glass state; there is 
essentially no choice. Since the glass curve moves to low 
density, this is a serious restriction. However, we may 
choose to work within the two-phase low-density fluid- 
crystal co-existence regime, but to the left of the glass 
curve (point marked B in Fig. ^). The outcome is then 
a question of kinetic control, and will not be completely 
settled by diagrams such as we are drawing. However we 
can make some educated comments. Thus, if we work 
in the two-phase regime of gas-solid {B), we may nu- 
cleate and grow crystals. Whether the proximity to a 
meta-stablc critical point is advantageous or not, as dis- 
cussed by Hagen and Frenkel is not our primary con- 
cern here, though this is an interesting proposition. The 
broader point is that by nucleating to the left of the glass 
curve one may enter the crystalline region (thus form a 
macroscopic crystal) by a route not described by the "adi- 
abatic" description here, and thereby avoid some of the 
complications of the glass. This renders the formation of 
crystal at least feasible, although where the glass curve 
runs through the two-phase region, it will remain difficult 
to form truly large high quality crystals. From this region 
where it is possible to crystallize, one should also exclude 
the two phase region, whether it is meta-stable or not, 
since the partial (micro-)phase separation, crystallization 
and glassification, all competing dynamically is unlikely 
to produce a good crystal also. This leaves only that re- 
gion bounded to the right by the glass curve, to the bot- 
tom by phase separation, and to the top by the gas side of 
the gas-crystal phase-coexistence as a likely candidate for 
forming good crystals. This is interesting. It leads us to 
suppose that for a fixed short range of the potential there 
is a "practical crystallization region" in the temperature- 
density plane, irrespective of the specific features of the 
equilibrium phase diagram. However, more importantly, 
there is also a limited regime of interaction ranges where 
such a slot is significantly large enough to be accessed 
experimentally. We believe that much of the discussion 
that has taken place in the last few years in the literature 
in relation to protein crystallization |]l|,p|j77| is almost cer- 
tainly heading in the correct direction. Thus, Mushol and 
Rosenberg [|j show in their Figure 12 the typical situa- 
tion for a globular protein phase diagram. They exhibit a 
phase diagram that has (a) two metastable liquid phases 
in equilibrium, a more and less dilute phase of protein 
(gas- liquid in our language), (b) a 'gellation curve' that 
we associate with the glass curves in our work, and (c) 
a fluid-crystal coexistence regime. They name the 'good' 
regime for crystallization zone I, and the others II and III. 
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We have also essentially partitioned our phase-diagram 
into the same types of zones, and concluded that this 
gas-crystal region described above by the glass, binodal 
and gas side of the coexistence boundary (zone I in their 
language) would be the most favourable for formation of 
crystals. For the sake of comparison in one of our figures 
(Fig.^) that has a range typical of globular proteins, we 
have marked regions I-III in analogy with rcf . ||^ . 
We would argue that the present work, with its phase 
diagrams, and discussion of control parameters is the 
quantitative expression of these ideas that have arisen 
in the protein crystal literature. This is potentially en- 
couraging, since it opens the possibility to make more 
quantitative study of these systems. 
However, perhaps the most promising directions involve 
the study of the current model, and underlying ideas in 
more depth to see what independent kinetic routes exist 
to the formation of crystals. We have seen how these 
phase diagrams are indicative of the kinetic behaviour, 
but we believe that there will be much more significant 
insights as to how to form high quality crystals in these 
regimes if a deeper understanding of kinetics is acquired 
in future. In particular the glass analogy seems promis- 
ing as a means to characterize the more confined phase 
space experienced by these systems. Aging, and kinetic 
phenomena in general, is an arena that is growing in im- 
portance |78-82| and may offer significant advances. The 
traditional viewpoint of activated processes, and sim- 
ple kinetic processes is without doubt incomplete in the 
limit where we approach the rather confined phase-spaces 
characterized by approach to a glass transition. The re- 
alization that we are in a 'glassy' scenario may well assist 
in the development of new theories of kinetics of crystal- 
lization more appropriate for such questions. 
In any case, one can hardly doubt the high degree of prac- 
tical significance that kinetic phenomena associated with 
short-ranged potential systems will have in the coming 
few years. Given that we discuss a model potential that 
is only slightly different from those long considered in 
liquid-state theory, we must be intrigued by the novelty 
in supposedly simple situations. 
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FIG. 1. Liquid-liquid phase diagram from SCOZA com- 
pared to simulation data for b=1.8. The GEMC data are 
taken from Shukia |Q (see text for details). The estimation 
of the critical points are from generalized means spherical 
approximation (GMSA), modified hypernetted chain (HNC) 
and HRT 
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FIG. 2. Liquid-liquid phase diagram from SCOZA com- 
pared to simulation data for b=6. The GEMC data are taken 
from Shukia for A'^ = 216 pq] (see text for details). 




0.3 0.4 

FIG. 3. Phase diagram for the Yukawa fluid with screening 
parameter b = 5.0. The crosses represent fluid-solid phase 
transition, the continuous line is the binodal and the dashed 
one is the spinodal. The filled circle is the critical point. The 
glass transition line as evaluated for Mode Coupling Theory is 
also displayed (open circles). The glass line shifted to obtain 
asymptotic value for T — » oo to be the experimental packing 
fraction (ji = 0.58 is presented (stars). The subscripts c and p 
refer to the critical point and the triple point respectively 




FIG. 4. As in Fig. | with b = 6.05. The fluid branch of the 
fluid-crystal coexistence line now passes trough the liquid-gas 
critical point. The shifted glass line is not represented in this 
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FIG. 5. As in Fig.| with 6 = 7.5. 
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FIG. 6. As in Fig.g with h = 30. For the points A and 
B see the text. At this h value the glass line passes trough 
the metastable liquid-gas critical point. The labels I, II and 
III are chosen by analogy with the proposition of Mushol and 
Rosenberg See text for details. 
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FIG. 7. As in Fig.g with 6 = 60. The subscripts c' and p' 
refer now to the the critical and triple points of the solid solid 
transition. The shifted glass line is not represented in this and 
in the following cases. 



b=100 



0.8 



T , 



0.6 « — 
H - P' 

CO 

0.4- 



0.2- 

T*- 

c 

0, 







X ^x 
XT 



0.2 



0.4 



0.6 



FIG. 8. As in Fig.H figures for h = 100. 
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FIG. 9. Phase Diagram for the SW model for A = 0.03. 
The crosses represent the sohd-fluid phase coexistence and 
the set of open circle is the glass line. Note the solid-solid 
coexistence on the high densities side of the phase diagram: 
its critical point is labelled by a filled circle. The position of 
the liquid-solid-solid triple point is also displayed {(j)pi,Tpi) 
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FIG. 10. The critical temperature Tc plotted as func- 
tion of the screening parameter b. The glass transition and 

the solid-fluid coexistence temperature at the critical packing 
fraction, 4>c, are also displayed. For completeness in the inset 
the d>c as a function of b is also shown. 
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